Manipulating motional states by selective vibronic interaction in two trapped ions 
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We present a selective vibronic interaction for manipulating motional states in two trapped ions, 
acting resonantly on a previously chosen vibronic subspace and dispersively on all others. This is 
done respecting technical limitations on ionic laser individual addressing. We discuss the generation 
of Fock states and entanglement in the ionic collective motional degrees of freedom, among other 
applications. 



PACS number(s): 42.50.Vk, 32.80. Pj, 42.50.Dv 

Coherent manipulation of vibronic states of two or 
more trapped ions has become a subject of increasing in- 
terest in the last few years [jD . A number of experimental 
and theoretical advances has been done recently in this 
field, aiming to make applications possible and to test 
basic features of quantum mechanics. In this context, 
Bell states and multiparticle entanglement of the ionic 
internal degrees of freedom have been studied and 
partially realized |^|j6). Beside their fundamental impor- 
tance, such states are essential ingredients for the imple- 
mentation of several interesting applications, as quantum 
logic and quantum computing using trapped ions |l],[7]|§ • 
One of the main problems in realizing some of the above 
proposals is the technical difficulty in the individual ionic 
laser illumination. Whereas many efforts are being done 
to overcome this problem H|l(| , it still represents an ob- 
stacle to the existence of realistic proposals for the deter- 
ministic manipulation of the motional degrees of freedom 
of two or more trapped ions. 

Recently, M0lmer and S0rensen |2) have shown that if 
one illuminates two two-level ions with two lasers that 
act dispersively in the blue and red first side band tran- 
sitions, it is possible to find an effective interaction, in 
the Lamb-Dicke limit, that causes a direct transition be- 
tween the ground and the doubly excited electronic levels 
of the two ion system, without changing their vibrational 
state. This permits the deterministic generation of en- 
tangled states that are linear combinations of these two 
electronic states. In ref. P], the combination of this inter- 
action with a carrier transition permits the generation of 
all electronic Bell states. This leads to more general ap- 
plications, as teleportation and entanglement swapping 
jLjJ . All these proposals are mostly concerned with the 
deterministic manipulation of internal states of the ions. 
However, it is also important to find ways for manipu- 
lating the vibrational degrees of freedom of two or more 
trapped ions, respecting the limitations imposed on their 
individual addressing. 

In this article we present an effective interaction con- 
sisting of two dispersive Raman pulses simultaneously 
illuminating two trapped ions, which opens the possibil- 
ity to manipulate coherently their vibrational motion. As 
the electronic Stark shifts, induced by the Raman pulses, 
depend on the motional state of the ions, the resulting dy- 



namics can be described by a Jaynes-Cummings-like in- 
teraction acting distinctly on different subspaces of their 
vibronic Hilbert space. As we will show below, the fre- 
quencies of the two dispersive Raman pulses may be cho- 
sen in such a way that the effective interaction becomes 
resonant to a previously chosen vibronic subspace while 
remaining non-resonant to others. This enables us to 
excite selectively a desired subspace inside the motional 
Hilbert space of the ions. Vibrational state engineering 
of both center of mass and relative motion can then be 
done by means of this special property. 

We consider two two-level ions of mass m, confined 
to move in the z direction in a Paul trap. They are 
cooled down to very low temperatures [ p"2[ and may 
perform small oscillations around their equilibrium po- 
sitions, zio = —d/2, Z2Q = d/2. We denote by Z = 
(zi + Z2)/2 and z = [z\ — ^)/2 the center of mass 
and relative position operators, respectively. Both ions 
are simultaneously illuminated by two classical homoge- 
neous Raman effective pulses Ej = E ie l ( q ' 1 ' r ~ UJlt ' > and 
En = Equ e t ^ 2 ' r ~ lA)Ilt \ with wave vectors q\ = q_2 = q, 
parallel to the z direction and angular frequencies loi and 
mn. The Raman pulses frequencies will be chosen to be 
quasi-resonant with a long-living electronic transition be- 
tween two ionic hyperfine levels | |j) and | lj) (j=l,2), 
with energy Hloq and 0, respectively. The total Hamilto- 
nian of the system may be written, in the optical RWA 
approximation, as 

H = H + H inU (1) 

with 

Hq = hu> (S-\-iS-i + S+2S-2) + hua'a + hv r Pb, 
H int = HSl e lq z (s +1 e iq& ' 2 + 5 +2 e - 4 « 2 / 2 ) x 

(e"^ 1 ' + e -^ Ijt ) + H.c. (2) 

Here S+j = |tj)(lj| is the flip operator associated with 
the electronic transition | lj) — > \ lj) in the ion j. The op- 
erators a and b (d^ and b ') are the annihilation (creation) 
operators associated with the center of mass mode of fre- 
quency v and the relative vibrational mode of frequency 
v r , respectively. For simplicity, we have assumed that 
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the same Rabi frequency fl (taken as real) is associated 
with both lasers. 

We start by taking the frequencies u>j and ion as: 

u>i = ujq + kv + k r v r — 8 
ujii = u Q + S, (3) 

with S <C v, v r . The frequencies are chosen so that ujj + 
u>n = k v + k r v r + 2co>o- 

Following the usual treatment for one single ion in- 
teracting with a laser field p3| , we make the rotating- 
wave-approximation (RWA) with respect to the vibra- 
tional frequencies and select the terms that oscillate with 
minimum frequency. In the interaction picture we obtain, 
for the interaction Hamiltonian, the following expression: 



6L = hn [S'a) k P K H kk e M + SlH ne 
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where 



Here 



S' + = S+ie^"/ 2 + (-) fe "5 +2 e- 
S'l = S +1 e l4>o/2 + 5 +2 e-^ o/2 
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= Zod is the phase difference due to the equilib- 
rium distance between the two ions, which, in this paper 
may be taken as 2mr without loss of generality , and 
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The parameters rj = q^h/Am v and r\ T = qy/Ti/Am v r 
are the Lamb-Dicke parameters corresponding to the CM 
and the relative modes, respectively. L™ are associated 
Laguerre polynomials, whereas the states \n) and \n r ) 
are the eigenvectors of the number operators a^a and 
b^b, respectively. 

If S is large enough so that S 3> f2 and for times t 
such that 5t ^> 1, it is possible to derive a two-photon 
effective time-independent Hamiltonian, where both CM 
and relative modes can be excited. It can be written as 
a sum of three terms: 



H e fi — H i + Hi + H 3 . 

with 
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S^S +j (a k b k -aH^Hl kr H^ Q ) + H.c. , (10) 



where Hq = fl 2 /5. The first term, Hi, gives rise to an 
anti Jaynes-Cummings dynamics, leading to a simultane- 
ous excitation (or de-excitation) of the electronic states 
of the two ions, accompanied by the creation (or anni- 
hilation) of k vibrational quanta in the CM mode and 
k r vibrational quanta in the relative mode. The factor 
e kr = (1 + (— ) kr ) prevents excitations of odd number 
of quanta in the relative mode, so that the symmetry 
by exchange of the two ions is maintained. The second 
term, H 2 , generates a dynamics where, simultaneously, 
one ion undergoes a transition from the ground to the 
excited electronic state and the other ion makes a tran- 
sition in the inverse direction. This process is not ac- 
companied by any excitation of the vibrational modes. 
The third term, H3, generates motional dependent dy- 
namical energy shifts in the electronic levels. Due to 
this dependence, this term turns the processes induced 
by Hi and Hi more or less resonant, depending on the 
particular level of excitation of the vibrational modes. 
Note that the sensitivity of the energy shifts to the vi- 
brational state of the ions increases with increasing val- 
ues of the Lamb-Dicke parameters. For not too small 
values of the these parameters, it is possible to make the 
interaction Hi completely resonant inside a previously 
chosen subspace {| j, I, N, N r ), \ T, j, N + fc, N r + k r )} , 
whereas remaining largely non resonant inside other sub- 
spaces. The same applies for Hi in a given subspace 
{| I, t, M, M r ), I t, j, M, M r )}. This may occur, for ex- 
ample, when the Raman lasers frequencies, originally 
given by Eq. ^, are modified and correctly tuned to take 
into account the motional dependent energy shifts. 

To describe these effects in more detail, we turn our 
attention to specific cases. We first consider excitations 
to the first blue side band of the center of mass mode 
(k = l,k r = 0). If we start from the electronic ground 
state, only Hi and H3 will be effective , and we may 
write: 

H cS = nn Gl{2{iri)S +1 S +2 {a}Fi,F ] + 
^(S +1 S_i + S + iS-i){r, 2 a)F 2 a - F 2 ) 

— (S-xS+i + S-iS + i){rt 2 aa)F 2 - F 2 ) + H.c.} . (11) 

As can be easily seen from Eq.|ll], the energy shifts of 
levels j I, I, n, n r ) and | |, |, n + 1, n r ) are given by: 

A?f ' = 2hn g 2 (n r ) [f 2 (n) - V 2 (n + l)f 2 (n)] (12) 
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respectively. In Eq. (|12| 



Aj+ lin ' = 27iQo<?oK) Win + l)f 2 (n) - f 2 (n + 1)] 



f k {m) = e 
9k (m) = e-tf/a 
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By properly adjusting the laser frequencies we may put 
them in resonance with the Stark shifted levels associated 
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with a previously chosen vibronic subspace {| J., J., N, N r ) 
, | f , |, N + 1, -/V r )}, while preventing resonant transitions 
in other subspaces with n N and n r ^ N r . This can 
be done by modifying the laser frequencies to 

ujj = & + v - 5, cjii = &q + 5, (14) 

where 2ujq = 2luq + A^ +1 ' Nr — A^ Nr is the renormalized 
splitting of the levels. Notice that, for very small values 
of the Lamb-Dicke parameters, the motional dependence 
of the dynamical Stark shift disappears. For this reason 
it is important to work beyond the Lamb-Dicke regime 
in order to effectively select a chosen subspace out of 
the whole vibronic Hilbert space. It is noteworthy to 
mention that for special values of the Lamb-Dicke pa- 
rameter 77, it may happen that energy shifts A™^~ ' nr and 
A"j™ r are equal for certain values of n, irrespectively of 
the state of the relative vibrational mode. For example, 
for 77 ~ 0.51, 0.42, 0.24, transitions inside the subspaces 
{| N,N r ), I i,1,N + l,N r )}, become resonant for 
N = 1,2,8, respectively. Clearly, in this case it is not 
necessary to correct the laser frequencies. 

In order to check this model, numerical simulations 
were done using the time dependent Hamiltonian given 
in Eq. (|). Starting from the state | |, 1,(^0) Jl| and 
selecting the laser frequencies as in Eq. (|14|), we were 
able to observe complete Rabi oscillations between the 
states I |,|,0,0) and | 1,1,1,0), in agreement with the 
model discussed above (see Fig. 1). In particular, for a 
tt pulse, the state | f,!, 1,0) is generated with « 100% 
efficiency. For the same laser frequencies, we plot the 
Rabi oscillations of the population corresponding to the 
initial state | |, j, 1, 0). We can see that, indeed, the Rabi 
oscillations for this case have a very small amplitude. 

Similar results may be obtained also for transitions 
leading to even excitations of the relative vibrational 
mode only. This is done by choosing k — in Eq. [)] and 
correctly tuning the lasers to take into account the self 
energy terms. We checked numerically that it is possible, 
in this case, to drive resonantly Rabi oscillations inside 
selected subspaces {| j, |, N, N r ), | f , T> N, N r + k r )}. 

The Hamiltonian (|9|) could be used to generate a large 
set of motional states. For example, any Fock state as- 
sociated with the center of mass motion and with even 
excitations of the relative vibrational mode could be ob- 
tained from the initial state | |,|,0,0) by successively 
applying tt pulses with different frequencies. However, 
if one is interested in generating a highly excited Fock 
state, this process could take an unsatisfactory long time. 
A more efficient, non unitary, way of producing such 
states, as well as engineering other vibrational states, is 
to start from a product of the electronic ground state and 
any motional state IV'vib)- By selecting the laser frequen- 
cies, we excite only a chosen vibronic transition {fc, k r } 
with a 7r pulse. Ideally, we would end up with a su- 
perposition of the two states | f, |, N + k, N r + k r ) and 
1,1) ® (IVVib) - |iV,iV r XiV,iV r |^ vi b)). Measurement of 
the electronic levels projects out either the Fock state 



\N + k, N r + k r ) or the original state with a "hole" in the 
N, N r component. 

As a numerical example we have taken the initial state 
to be the product of the electronic ground state with a 
coherent state of the CM motion and the vacuum of the 
relative motion. We start from a coherent state |lj] with 
n = 4.0 and excite the center of mass transition from n = 
4 to n = 5. A dark event in a fluorescent measurement 
should leave us with an state that is close to the Fock 
state with n = 5. In Fig. 2 we show the results for the 
phonon distribution. As expected, it is possible to make 
a "hole" in the vibrational quanta distribution ( Fig. 2a 
) while creating a quasi Fock state ( Fig. 2b) using a tt 
pulse. The state with a "hole" will be associated with the 
ground electronic state, while the approximate Fock state 
of Fig. 2b will be associated with the excited electronic 
state. As shown in Fig. 2b, small contamination occurs 
around the target Fock state because transitions to levels 
other than n = 5 are not totally suppressed. For the case 
studied, the probability of finding the approximated Fock 
state after fluorescence is about 30%. A Fock state with 
n = 6, for example, can now be obtained if we apply 
subsequently to the ions another tt pulse resonant to the 
transition | |,f,5,0) <-»-| 1, 1, 6,0). If by measuring the 
electronic states we find | J., j) ( a priori probability of 
« 87%), we are left with a state very close to the Fock 
state |6, 0). In this case the fidelity is «s 99% (See Fig. 2c). 
Depopulation of a region around a certain value of N, 
can also be achieved by proper choice of the Lamb Dickc 
parameter and the detuning 5, taking advantage of the 
quasi resonance character of the interaction for n's very 
close to N. 

Another interesting application is the generation 
of maximally entangled states inside the subspace 
1,0), |0,2), |1,2), |0,0) of the vibrational modes. Start- 
ing from the electronic and vibrational ground state, a 
7r/2 pulse with a given frequency would generate the 
state [| i, I, 0, 0) + I t, T, 0, 2)]. A subsequent tt pulse 
with another frequency would then lead to the state 
-^\ f,T) [| 1 , 0) + |0, 2)] which is a maximally entangled 
state. Similar procedures could also lead to any maxi- 
mally entangled state. 

Entanglement transfer may also be achieved between 
the internal and motional degrees of freedom 
For example, assume that we start with the state 
-73 [| I) + I T, T)] ® |0, 0). If we apply to the ions a tt 
pulse, connecting the states | J., J. 0, 0) — > | f, 1,1,2), we 
obtain | f, |) eg) (|0, 0) + |1, 2)). In this case the efficiency 
is very high since the state | \, "f ) <g> |0, 0) does not couple 
to any state in an anti Jaynes-Cummings transition. Of 
course, the inverse process, where entanglement is trans- 
fered from the motional degrees of freedom to electronic 
one, is also possible. 

In summary, we have engineered an interaction that, 
respecting limitations of ionic individual addressing, en- 
hances our possibilities of manipulating and generating 
diverse vibronic states of two trapped ions. This interac- 
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tion acts selectively in a previously chosen vibronic sub- 
space, {| l,lN,N r },\ ttN + k,N r + k r )}, permitting, 
in principle, a complete transfer of populations inside this 
subspace. The applications mentioned above are only a 
few relevant examples of what may be done by selectively 
addressing the initial vibronic states. 
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FIG. 1. Center of mass first side band blue excitation. 
The larger oscillations correspond to the population of the 
state | |,|,0,0). The smaller oscillations correspond to the 
population of the state | j, J., 1,0). rj — 0.5, S — 40r]£l. 



FIG. 2. Figure 2. Vibrational distributions P(n) for the 
center of mass motion after a 7r pulse resonant to the tran- 
sition from n = 4 to n = 5. The initial state is the product 
of a coherent state with n = 4.0 and | j, J.), (a): Distribution 
correlated to | |,|). (b): Distribution correlated to | T,T). 
(c) Distribution correlated to | J., J.) after the second pulse. 
f] = 0.3 and S = 40.0r?fl 
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